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29 . ABSTRACT 

o 
> 

^ ' Some years ago, a new powerful technique, known as the Classical Effective Field Theory, was 



proposed to describe classical phenomena in gravitational systems. Here we show how this ap- 
proach can be useful to investigate theoretically important issues, such as gravitational radiation 
in any spacetime dimension. In particular, we derive for the first time the Einstein-Infeld-Hoffman 
Lagrangian and we compute Einstein's quadrupole formula for any number of flat spacetime di- 
mensions. 
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B Useful relations for Feynman diagram computations [28 



1 Introduction 

One of the most conceptually simple and interesting problems in any theory of gravitation concerns 
the dynamics of two structureless particles, moving under their mutual gravitational interaction. 
This is the two body problem [T], which has proved to be fairly simple in Newtonian gravity but 
a formidable task in General Relativity. A clear understanding of its solution allows for tests of 
Einstein's and alternative theories of gravity [21 [3], in the weak and strong field regime. With 
the advent of high-sensitivity gravitational wave experiments, such as LIGO and others [1], the 
importance of the two body problem has exponentially increased. Indeed, one of the main targets 
for these experiments is the inspiral and coalescence of neutron star or black hole binaries. An 
accurate solution to this two body problem is necessary not only to compare theory with experiment 
but also to improve the chances of detection: the signals expected to impinge on the detectors will 
be completely buried in noise, and the most powerful way to dig them out, matched filtering 
[5], requires a very accurate knowledge of theoretical templates. This curious interplay between 



theoretical and experimental needs has boosted the effort to find accurate solutions to the dynamics 
of gravitationally interacting, but otherwise isolated bodies in General Relativity. 

Only recently have long-term stable numerical evolutions of black hole binaries been achieved 
[6], but these still cover only a very restricted part of the evolution, since they are extremely 
demanding from a computational point of view. Even if the situation improves dramatically in 
the near future, there's only a limited amount of computational power and therefore numerical 
waveforms always cover only a limited amount of time in the dynamics of the two body problem. 
Analytical solutions are required, and they serve two very important purposes. First of all, analyt- 
ical solutions are valuable on their own, as solutions which in some regime are accurate enough and 
describe well enough the problem at hand. They usually provide physical insight that is absent on 
numerical calculations. A second important purpose relates to the validation of numerical codes, 
allowing simultaneously a extension of the numerics by some kind of matching procedure. Ana- 
lytical solutions to the two body problem have been studied for many decades now, starting from 
Einstein himself. The traditionally most powerful technique, the Post-Newtonian (PN) formalism 
[3 El E], expands the problem as a series in powers of vjc. 

In an impressive tour de force, Goldberger and Rothstein [lOl [TTl [12] and Porto and Rothstein 
|13j-|18j and later Kol and Smolkin [19] have shown how to extend and simplify considerably the PN 
calculations, using methods borrowed from the Effective Field Theory (EFT) approach to quantum 
field theory. As discussed in [10l[19], this EFT of gravity is a classical theory, justifying the use of 
the acronym CIEFT (Classical Effective Field Theory) [19]. Although the computations are carried 
out in terms of Feynman diagrams using the quantum field theory language, no quantum effects 
are considered. The objects whose physical properties one is interested in are macroscopic, and 
hence the quantum corrections to the classical observables are negligible. This method is especially 
well suited to treat the early stages of the binary inspiral problem when the two objects are non- 
relativistic. In this problem there are three widely separated scales: the typical internal scale of 
the objects tq, the orbital distance r between them, and the radiation wavelength i ~ r/v, where 
V is the typical velocity of the system. At the computational level, having an exact description 
of the system does not seem to be practical nor necessary. The reason is that in gravitational 
wave physics one is only interested in the radiation that reaches the detector, and therefore it is 
sufficient to have an effective description of the system that captures this information. This is 
precisely what the CIEFT does: by systematically integrating out short distance scales, one is left 
with an effective description of the system in terms of the relevant degrees of freedom at the scale we 
are interested in. For the binary problem, these are radiation gravitons coupled to point-particles, 
and the relevant length scale is i. The main advantage of this approach is that these degrees of 
freedom are described by a Lagrangian and to obtain the observables one has to compute Feynman 
diagrams, as in standard quantum field theory. Therefore, the procedure is very systematic. 

Another advantage of the CIEFT approach is that it has manifest power counting in the small 
parameter of the theory, namely, v/c. Therefore, one can systematically compute a given observable 
to any desired accuracy since one can determine which terms in the Lagrangian and which Feynman 
diagrams will contribute to that order. This is important from a practical point of view because 
highly accurate theoretical predictions are required to compare to the putative gravitational wave 
signals to be detected by the interferometers. The CIEFT approach to gravitational radiation also 
provides a systematic way to deal with the divergences that plague the PN calculations. These 
divergences can be thought of UV divergences that arise because one is treating extended objects. 



like black holes or neutron stars, as being point-like. In EFT these divergences can be simply dealt 
with dimensional regularization. Moreover, the finite-size effects can be easily taken into account 
by including in the Lagrangian higher order operators that are consistent with the symmetries of 
the problem. The coefficients of these operators are then determined by a matching calculation in 
the full theory, see ^0\ |T2] for a more detailed discussion. 

In most realistic situations one expects that the two compact objects in the binary system are 
spinning. Therefore, from both a theoretical and experimental point of view it is important to 
incorporate rotation to the CIEFT. This has been done in |13|-|2U|. The resulting EFT of gravity 
has manifest power counting and one can compute observables using Feynman diagrams. Therefore, 
it retains the computational power of the original CIEFT of gravity of [lOj . 

Finally we should point out that the CIEFT approach has also been successfully applied to 
other problems of interest which involve widely separated length scales. In particular, it has been 
shown that the CIEFT methods are particularly useful to compute, in a perturbative expansion, 
the physical parameters of caged black holes |19tl21j. In this situation, there are two length scales: 
the size of the black hole rg and the length of the compact circle L, and the small parameter is a 
power of their ratio. 

What we want to show here is that the CIEFT approach can and should be used in other 
contexts as well, in particular in calculations of gravitational wave emission in higher dimensional 
scenarios. Gravitational radiation properties (intensity, distribution, etc) depend on the dimen- 
sionality, curvature, etc., of the model under consideration. Thus, the study of gravitational waves 
in different higher dimensional scenarios may provide a means to constraint or distinguish different 
models. Wave propagation in higher dimensions is intrinsically different from the four-dimensional 
case (221 [231 IMl ISSl EH 123 EH] : waves in flat even-dimensional spacetimes propagate on the light 
cone, just as in our usual four-dimensional world, but waves in odd dimensional spacetime do 
not. The extension of the notion and properties of gravitational waves to any even-dimensional 
spacetimes is almost straightforward. By isolating the most important contribution to the Green's 
function, we recently computed the field in the wave zone and the quadrupole formula in general 
even-dimensional spacetimes ^27j. This formula provides the simplest, lowest-order contribution to 
gravitational radiation emission in general situations of physical interest. It is relevant, for instance, 
for string-motivated braneworld scenarios where string-field gravity plays an important role |31j . 

Unfortunately, wave propagation in general odd-dimensional spacetimes are much harder to deal 
with [221 \29\ I30j . In particular, some of the rules we are used to are not valid in these spacetimes. 
For instance, in odd-dimensional fiat spacetimes a propagating wave leaves a "tail" behind. This 
means that a pulse of gravitational waves (or any other massless field) travels not only along the 
light cone but also spreads out behind it, and slowly dies off in tails. Progress in the study of 
gravitational waves has been hampered by this fact. 

We will show in this work that the CIEFT approach is able to give us the correct quadrupole 
formula in fiat, even dimensional spacetimes and in the same stroke extend that formula to odd- 
dimensional spacetimes. The fundamental reason for this, other than the simplicity of the CIEFT 
approach is that it works in the momentum space, where the Green's function have the same 
functional behavior. At the same time, we will also derive for the first time the generalization of 
the Einstein- Infeld-Hoffman [32] Lagrangian to an arbitrary number of spacetime dimensions. 

The plan of the paper is the following. In section [2] we start by describing the problem and 
review the CIEFT approach [10] that we use. In section [3] we find the higher dimensional Einstein- 



Infeld-Hoffman Lagrangian. Finally, in section [J] we determine the quadrupole formula in a general 
flat higher dimensional spacetime. In appendix |Al we derive the Feynman rules for the graviton 
propagators, graviton vertices and point particle vertices needed in the main body of the text. In 
appendix [B] we present some useful relations needed to compute Feynman diagram contributions. 
We try to be self-contained in our presentation. 

We take d to be the number of spacetime dimensions. Newton's constant in d-dimensions is 
represented by Gd and is related to the gravitational coupling Kg through (13. 4p . We take the 
(+, —,•••,—) signature for the Minkowski spacetime. The speed of light is taken to be c = 1. 

2 Description of the problem. Effective field theory approach 

2.1 Lengthscales in a binary system 

We are interested on finding the first PN correction to the gravitational interaction between two 
bodies, and the quadrupole formula expressing the gravitational energy emitted by a system moving 
at low velocities. One of the simplest background to address simultaneously these issues is a binary 
system. The interacting bodies can be black holes or neutron stars with typical mass m. Such 
system has three lengthscales. To start with, we have the size rg oc m of the gravitational bodies. 
Then, we have the typical length separation r between them. Finally, we have the wavelength i of 
the emitted radiation. Alternatively, we can make use of the typical velocity v of the system. Then, 
its orbital period is T ~ r/v, and its evolution sources the emission of gravitational waves with 
frequency T~^ and wavelength i ~ r/v (here and henceforth we take unit light velocity, c = 1). If 
the system's constituents move slowly compared with the speed of light, u <C 1, we clearly have 
three widely separated lengthscales, 

ro^r^i. (2.1) 

In this slow motion condition we can study some important properties of the system using a post- 
Newtonian (PN) expansion in powers of w <^ 1. Alternatively, as we will do, we can follow the 
CIEFT approach, specially tailored for these problems. Both formalisms consist on systematically 
solving the Einstein's equations with non-relativistic (NR) sources by taking a power series expan- 
sion in the small velocity parameter. The three lengthscales are not independent. Indeed the Virial 
theorem for a small motion system governed by an almost Newtonian interaction states that 
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-^ ~ ;f3 ' (2-2) 

where k^ is the gravitational coupling proportional to Newton's constant Gd', see (13. Ih . This relation 
between the kinetic energy and the Newtonian potential will be important later since it tells us 

2 

that -§:^ also contributes at the same order as v"^ in a velocity expansion. 

The gravitational field gjj,u created by this system is a solution of Einstein's equations, which 
follow from extremization of the Einstein-Hilbert action 

SEH[gf] = \ [d^'x^Rigf], (2.3) 

where ^/gj is the determinant of the metric gfn, and R[gf] the corresponding Ricci scalar. As is well- 
known, this yields a non-trivial non-linear system of differential equations and no exact solution 
gfiu is known for the two-body problem. 



2.2 Integrate out the internal structure or short lengthscale physics 

In the CIEFT framework we first decompose the full, or exact gravitational field gfi^ into two 
components. Take g* to be the metric that describes the gravitational field in the short lengthscale 
region of order rg, and Qf^i, the metric encoding a similar information but this time in the region 
with lengthscales of order r and larger. Then, within a good approximation, we can decompose the 
full metric as 

dU = alu + ffM^ • (2.4) 

We can now do a first disentanglement of lengthscales, namely we separate the short lengthscale 
gt^jj from the long lengthscale g^jy. Within a traditional effective field approach, this is achieved by 
integrating out 5^^, through the path integral 

exp {iSeff[g^,u]/h) = I Vg'^^ exp (iSEH[9l,.,x]/hj , (2.5) 

where h is Planck's constant. In the classical limit, /i — > 0, this path integral reduces in the saddle 
point approximation to the computation around the classical solution. This is the regime of interest 
for us here, since we want to use EFT techniques to compute purely classical (tree-level) results 
(for a clear discussion of this issue see [19] ) . In this classical limit we can then write 



where 



Seff[g, Xa] = SEH[g] + SmIq, Xa], (2.6) 

SEH[g] = [ d^x L[g] , L[g] = \ ^gR\g\ ; 

ShAg, Xa\ = Spp[g, Xa]-\ , Spp[g, Xa] = -'^rua / dTa . (2.7) 

a '' 

This effective action describes the motion of the point particle ensemble in the g^i, background. 
It preserves the symmetries of the original theory, namely general coordinate invariance, world- 
line reparametrization invariance and, since we do not consider internal spin degrees of freedom, 
SO{d — 1) invariance. The black hole or neutron star effective action Sm is to leading order 
given by the point particle effective action Spp, where dTa = ^J gp.v{xa)dxadx'^ is the proper time 
along the worldline Xa of the a**^ particle. The dots in Sm represent higher order terms that 
describe non-minimal couplings of the point particles to the spacetime metric and are discussed 
in detail in [12]. They account for finite size effects. Therefore they are not relevant to com- 
pute the Einstein-Hoffmann-Infeld correction to the Newtonian interaction neither to determine 
the quadrupole formula we are interested in. Indeed, these quantities are clearly independent of 
the internal or finite size structure of the gravitational bodies. They would be essential to address 
dissipative or absorption processes on the black hole horizon or star's surface |12[ I17j. 

Resuming, after the first disentanglement, where we integrated out the internal structure of the 
gravitational bodies, we have a low energy effective theory of point particles coupled to gravity 
described by (|2.6p . 



2.3 Potential and radiation gravitons 

We can now proceed to further disentangle the two other lengthscales of the problem. Since we 
are in the regime (j2.ip . i.e., r,i ^ vq we can treat the long lengthscale metric Qf^^, perturbatively. 
That is, for distances much larger than tq we are in the weak field regime and we can linearize Qfj^i, 
around flat spacetime rj^^, 

= Vfj-u + Kg {Hfj_u + hfj_u) , \Hf_tu\,\hfj,u\ <.l- (2.8) 



We have decomposed the small perturbations Sg^u into two components [10]. The so-called potential 
graviton if^,y describes the gravitational field for lengthscales of order of the orbital length r, while 
hf^iy is often denoted as the radiation graviton describing the gravitational field at larger distances 
of the order of ^ ~ r/v. Therefore, the potential graviton Hf^i, is naturally the mediator of the 
gravitational interaction between the point particles, i.e., it is responsible for the force that binds 
the 2-body system (or n-body). On the other hand, the long wavelength radiation gravitons describe 
the gravitational waves emitted by the system during its time evolution. They are the ones that 
can propagate to infinity and be detected by an asymptotic observer. 

In the regime of p.Sp . a small v expansion of the point-particle action Spp yields 

Spp = 5^ma I'dxl (2.9) 

2 

where we defined the velocity vector v^ = ^^, and we kept only the low order terms in the 
expansion. Higher order terms not shown will not be used. This action describes the non-linear 
interactions between the point particles and the gravitational field. 

Let's have a closer look at the properties of the two kind of gravitons. The potential gravitons 
Hfj^^ have momentum that scales as k'^ = {k^,]<.) ~ (r' r)' Indeed, their energy is given by the 
source's frequency but, since they have lengthscale r, their spatial momentum is considerably 
larger, |k| ~ -. Therefore, they have spacelike momentum, k^^k^ < 0. They are off-shell and 
cannot contribute to propagating degrees of freedom to infinity. As an important consequence, 
derivatives of potential gravitons scale as daH'"' = {dtH'"' ,diH'"') ~ (f,i) H^'''. That is, in a 
small velocity expansion the spatial derivative is one order lower in v than the time derivative and 
gives the leading order contribution. To emphasize this order distinction between time and spatial 
derivatives it is useful to work with the Fourier transform of H^i,, 

H,A^) = I e^-'^H.J.O) . l^J^^ (2.10) 

A time derivative acts as dfH -^ J dtH\^ with dt ~ v/r, but a spatial derivative gets now replaced 
by a spatial momentum factor, diH -^ Jkji/k with |kj| ~ 1/r. The advantage of working with 
the Fourier transform is that the counting of powers of v is now much more explicit: a term with 
a time derivative J dtH\^ is immediately identified to be one order higher in v than a spatial term 
j \<.iHi^. This provides a practical "visual" advantage when power counting a long expression. 



We now turn to the radiation gravitons. Their momentum scales as k'^ = {k^,k.) ~ (^, ^), i.e., 
since they have lengthscale r/v their spatial momentum has the same scale as their energy. So, 
these gravitons are on-shell, A;^/c^ = 0, and propagate at the speed of light to infinity where they can 
be detected. Derivatives acting on a radiation graviton introduce a power of v/r, dah^y ~ ^^i/, 
i.e., time and spatial derivatives contribute equally in the power expansion of v, contrary to what 
happens with the potential gravitons. Not less important, the CIEFT relies on a small velocity 
expansion. So, when power counting powers of v we must do a multipole expansion of the radiation 
graviton /i^iy(x),ll| 

where (5x = x — X and X = (^ maX„) / ^ nia is the center of mass of the system (or any other 
reference point of the particle ensemble). That is, we must do a Taylor expansion of h^^^x) 
around the center of mass to consistently identify the leading and subleading contributions in the 
V expansion. Indeed note that 6x^ scales as 6x^ ~ r and the derivatives scale as dih^i, ~ v/r. 
Therefore each new term in the multipole expansion contributes with an extra power of v relative 
to the previous Taylor term. 

Summarizing, the ingredients to an EFT with a well-defined power counting scheme in powers 
of V are the decomposition of the gravitational field into the potential and radiation contributions, 
the expansion of the particle action, the Fourier transform of the potential graviton and finally the 
multipole expansion of the radiation graviton. 

2.4 Integrate out the intermediate orbital lengthscale 

To derive an EFT that has manifest velocity power counting rules, we now need to integrate out 
the orbital scale, i.e., the potential modes H^^, by computing the functional integral 

exp{iSNR[h,Xa]/h) = 'DH^^^^ex.p{iSeff[g,Xa]/h) , h^O, (2.12) 

where Seff[g, Xa] is defined in (j2.6p and the computation is done in the classical limit h ^ 0. In the 
process, the radiation graviton is treated as a slowly varying background field. After this operation 
we are left with a new effective action S]yR[h,Xa\ that describes the non-relativistic (NR) or Post- 
Newtonian approximation to General Relativity, in an expansion in powers of v. It encodes the 
information concerning the gravitational interaction between the point particles, mediated by the 
potential gravitons, but also the coupling between the emitting system and radiation gravitons. 

The explicit operation of integrating out the potential gravitons will be done in the next two 
sections using the appropriate Feynman diagrams and techniques. The action S]yR[h,Xa], to the 
desired order in v, is given by a sum over Feynman diagrams that must satisfy some rules. 
i) To start with, note that the worldline particles will be represented by a solid line. The point 
particles worldline's represent very short wavelengths and thus "infinitely" heavy fields that were 
integrated out to get Seff[g,Xa]. They thus have no associated propagator (roughly, their propaga- 
tor would be TjjT — 7 ~ —J since they are heavy fields) and are treated as background non-dynamical 
fields. So in the Feynman diagrams we have no sum over the momentum of the point particles. 



^Indices are lowered with the Minkowski metric 77^1/ 



This justifies why loops that are closed by the particles are not quantum loops but give instead a 

contribution to the tree- level result. It also justifies that we should only consider diagrams that 

keep connected when we remove the particle worldlines. 

ii) Next, in the Feynman diagrams, potential graviton propagators can appear only as internal lines 

but never as external lines. The reason being that potential gravitons have interaction range of the 

order of only the orbital distance r. Thus, they mediate the gravitational interaction between the 

particles but cannot propagate to the asymptotic region. 

iii) Finally, the Feynman diagrams can only contain external radiation graviton lines that propagate 

to infinity. These gravitons have long wavelength and do not contribute to the binding force between 

particles. So, there are no internal radiation graviton lines. 

Summarizing, integrating out the non-dynamical potential gravitons generates the gravitational 
interactions between the NR particles. In the functional integral (|2.12p . the effective action for 
these interactions arise from diagrams with no external radiation gravitons. This action will be 
computed in section [3J To order 0{v'^) it is the Einstein-Infeld-Hoffmann Lagrangian, Leih- On 
the other hand, the functional integral (|2.12|) also provides the coupling of radiation gravitons to 
the particle ensemble through the diagrams with external radiation gravitons. This action L^ad 
will be computed in section 14.11 

To compute the Feynman diagram contributions we need the Feynman rules for the graviton 
propagators, graviton vertices and point particle vertices. In appendix |X] we give these rules and 
the details of their computation. To find these Feynman rules we have to introduce the standard 
gauge fixing and ghost Lagrangian contributions. We will follow the background field method 
|331 [3^ that fixes the gauge in such a way that preserves the invariance under diffeomorphisms 
of the background metric. This gauge fixing scheme guarantees that the obtained action is gauge 
invariant. So the gauge fixing action satisfying these properties, given in Eq. ()A.4p . must be added 
to the action in the rhs of Eq. (j2.12p . There is also a ghost field contribution, Eq. (|A.4p . but this 
only enters in quantum loop corrections, not considered here. Thus, the ghost contribution will 
not be considered in the functional integral (j2.12p . 

2.5 Integrate out the long wavelength radiation scale 

Having obtained the Lagrangian -L^ad describing the interaction between radiation and NR particles 
from Snr-, we can finally integrate out the radiation gravitons h^y through the functional integral 

ex.p{iSQF[xa]/h) = 'Dhfj_yexp{iSrad[h,Xa]/h) , h^O. (2.13) 

Integrating out hf^^ means that we are left with an effective action SqfIxo] that depends only on 
the particle coordinates. Therefore, this integration arises from Feynman diagrams that have no ex- 
ternal radiation graviton lines. It yields the desired quadrupole formula. The explicit computation 
will be done in section [321 



2.6 Power counting scheme 

The graviton field was decomposed into the potential and radiation components, (j2.8p . Further- 
more, we did a small velocity expansion of the particle action (j2.9p . we took the Fourier transform 
(|2.10p of the potential graviton, and we did a multipole expansion (j2.1ip of the radiation field. 
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Table 1: NRGR power counting rules for time coordinate x^, spacetime derivatives 9^, potential 
spatial momentum k, potential graviton H\f_ ^, radiation graviton /i^j^y and the coupling Kgm. Note 
that after the Fourier transform treatment (j2.1Up of the potential gravitons, derivatives of any 
graviton introduce for sure a factor of v/r. When power counting we should have in mind the 
useful relation, [ng] 



These operations are key steps in the CIEFT formalism since they provide the arena for a clear 
identification of power counting rules in the small velocity expansion. Having it at hand we can 
uniquely assign powers of the velocity parameter v to any Feynman diagram. We therefore can 
organize systematically the Feynman diagrams in powers of v and compute their contribution to 
our observable up to the order we wish. Before computing in next sections the observables we are 
interested on, it is useful to list these power counting rules. 

Let us start with the potential graviton. Its propagator, here represented in short as {Hi^^^H^^^) 
is given by ()A.13p . Time scales as x" ~ r/v and potential graviton momentum goes as |k| ~ 1/r. 
Since a delta function scales as the inverse of its argument, (JA.ISP tells us that the propagator for 
f^k^^ goes as {H\^^ Hf^^^) ~ [(^(k)] [k^^] [5(x'^)] ~ ^d-\^2yj^ ^ ^d^^ Therefore, this fixes the scale 
of the potential graviton as H\^^^ ~ ^1/2^(^/2 _ 

Next, consider the radiation graviton whose propagator (haph^u) is written in (jA.Sp . Since 
the radiation graviton momentum scales as \k^\ ~ v/r, this propagator scales as (/la/j/i^i/) ~ 
[d'^k\ [A;~^] ~ (v/r) {v/r)~ ~ ^2-d^d-2 ^^^^ ^j^^g ^j^g radiation graviton has scale h^^, ~ (f/r)~2~. 

Finally, to find the scale of Kgm (where recall that k^ oc Gd) we introduce the orbital angular 
momentum, 



L 



mvr . 



(2.14) 



and we use the Virial relation (j2.2p to get k^ m? ~ mv'^r'^~^ = Lvr'^~^. Thus, Kgm ~ L^''^v^''^r^^ . 
These power counting rules are summarized in Table [TJ Furthermore, in Table [2] we collect the 
power counting rules for the graviton correlation functions and point particle vertices that are 
derived and presented in Appendix El 

Summarizing, an inspection of these power counting rules reveals that tree-level classical re- 
sults always follow from Feynman diagrams with a power in the angular momentum of L^ (for the 
diagrams contributing to the gravitational interaction between bodies) and L^' ^ (for diagrams de- 
scribing the interaction of the system with radiation). Quantum loop corrections would correspond 
to diagrams with extra powers of h/L <^ 1. So, the orbital angular momentum is the parameter 
that counts loops in this classical EFT. Keeping attached to tree-level diagrams we can get our 
classical observables up to the desired order in the velocity u <C 1, having always in mind that 
loops closed by a "heavy" particle worldline are not quantum. 

The unique power counting scheme providing a straightforward systematic analysis, and the 
clear representation of the computation provided by the Feynman diagrams are probably the 
strongest qualities of the gravitational classical EFT formalism. On the top of this, the concepts 
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Table 2: Power counting rules for the: radiation graviton propagator (ha/sh^y), potential graviton 
propagator (ifka^-f^q^^), 0{v'^) correction to the potential graviton propagator (-ffk^^-f^qp^)®; 3- 
potential graviton correlation function {Hi^^^ Hq^^Hp_^^), 3-radiation-potential graviton correlation 

function {hHi^^^^Hq^^), and for the several point-particle vertices V^l, displayed in (jA.36p . These 
graviton propagators or correlation functions and point particle vertices are derived and presented 
in Appendix lAl 



and techniques of regularization and renormalization are naturally incorporated and can be used 
to fix divergencies that appear at the classical level. For example we will later used dimensional 
regularization. 



3 Einstein-Inf eld- Hoffmann Lagrangian in higher dimensions 

In Newton's gravity theory, perturbations propagate at infinite velocity and only the mass sources 
the force between gravitational bodies. Einstein realized that not only the rest energy but also 
kinetic and graviton energies contribute to the gravitational interaction of a system, and that 
gravitons propagate not instantaneously but at the speed of light. Newton's gravitational force has 
therefore General Relativity corrections that can be organized systematically in a power expansion 
in the velocity of the interacting bodies. The leading f;-dependent correction of order 0{v'^) has first 
computed by Einstein, Infeld and Hoffmann in 1938 [S^ (in 4 dimensions) using a PN formalism and 
reproduced in [10] and later in [20] using the classical EFT approach. In this section we compute, 
for the first time, the Einstein-Infeld-Hoffmann correction to the gravitational interaction between 
two bodies in any dimension using the EFT approach introduced in [lOj. 

As discussed previously, the potential gravitons govern the intermediate lengthscale physics at 
distances of the order of the orbital length r. Therefore they mediate the gravitational interaction 
between the NR bodies. The radiation gravitons, being long wavelength fields when compared to 
the typical distance between the particles, do not mediate this interaction. We thus have to compute 
the functional integral (j2.12p keeping only Feynman diagram contributions with (internal) potential 
gravitons and no radiation gravitons. 

In subsection 13. II we first compute the 0{v^) Newton interaction and fix our choice for Newton's 
constant. Then, in subsection 13.21 we find the Einstein-Infeld-Hoffmann correction. 
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3.1 Newtonian Lagrangian 

The Feynman diagrams that contribute to leading 0{Lv^) order are those of Fig. [TJ The first 
diagram has a non-vanishing contribution given by 

Fig. [a = V^^^{H^^^{xl)H^^^{xt))Vll^ 

= ^ j^ / dx\dx^2^{x\ - xl) / j^e '"^ ^""^ ''^> 

^*d-2i6^^|xi(t)-X2(t)|'^-3' ^•'■'> 

where we used the Feynman rule (|A.36p for the point particle vertex V,'t'^ and (|A.13p for the potential 
graviton propagator {H^^^^Hq^^). We also used relations (jB.ip - fjB.Sp . In the end we did the 

identification t = x^. Note that this diagram scales as Lv^ since Vti-\{Hif^Hq)V/i\ ~ (L^'^?;'^) ~ Lv^ 
(see Table ED. 





v« 






T 




*■ 


—l 


>■ 



a) b) c) 

Figure 1: Feynman diagrams that give the leading Newtonian interaction. They are of order Lv^. Actually, 
the self-energy diagrams b) and c) give a vanishing contribution after doing dimensional regularization. Solid 
lines represent particle worldlines and dashed lines the potential graviton propagator (|A.13[) . Dots are point 
particle vertices. In these diagrams they are V^,^ in (|A.36[) . Note that loops closed by a "heavy" particle 
worldline are not quantum and contribute to tree-level results {i.e., we have no integration over the particle 
momentum). Quantum loop corrections would correspond to diagrams with extra powers of h/L <C 1, i.e., 
with graviton loops. (The reader interested in quantum corrections to the gravitational interaction between 
two bodies can see [371|3H1I3S] and references therein.) 

The two last self-energy diagrams in Fig. [1] are pure counterterm and have no physical effect. 
They renormalize the particle masses and formally vanish after doing dimensional regularization, 
as pointed out in [lOj. For example, diagram [lb) gives 

Fig.mb = (-1)'^ (^) .ImijdxUxlj^ j^ , (3.2) 

where k is the momentum circulating in the loop. The last integral in the expression above is 
divergent and it gives rise to a renormalization of the particle's mass. However, in dimensional 
regularization it can be formally set to zero by noting that it arises as the x — > limit of the 
integral 

e"^"-^ 7T3v; = —1^ ^^„. ^ \^\ ■ (3-3) 



(2^)'^-! {k^r (4vr)^ r(a) V4 

11 



Sending x ^^ before fixing d shows tliat this integral vanishes. An obvious equivalent argument 
also makes Fig. [Ih = 0. 

At this point we specify our normalization for the gravitational coupling Hg in terms of the 
d-dimensional Newton's constant Gd- We choose, as in d = 4, to work with the normalization 



K^ = 327rGd . (3.4) 

Defining the gravitational coupling in this way {i.e., in a d-independent form) has the advantage 
that the black hole entropy in d dimensions is always horizon area divided by 4Grf. For this reason 
this is the normalization most commonly chosen (see, e.g., [35j . |27j . |36jrl. With this choice, the 
0{v^) Newtonian potential energy reads 



d — 3 2r(-2-) 0^17111712 

d^ ^^ |xi(i)-X2(t)h 



^A' - 3 ^ d^ — mT\ „ u\\d-3 ■ ^'^■'^> 



3.2 Einstein-Infeld-Hoffmann Lagrangian 

The Feynman diagrams contributing to the next-to- leading order 0{Lv'^) are those of Fig. [2j Their 
sum gives the Einstein-Infeld-Hoffmann correction to the Newton interaction. 

► ^ ► 1 ► Y ► 1 ► • ► 



(l'^^2) / \ +(l^e^2) 

/ \ 



V' V" V" v" 

a) b) c) 



/ 



^ +(l-e^2) 



v» 
► T ► 

. i ► 



V" V" V" 

d) e) 

Figure 2: Feynman diagrams for the Einstein-Infeld-Hoffmann correction to the Newton interaction. They 
are of order Lv'^ . Only (internal) potential gravitons contribute as mediators of the interaction. 

Using Feynman rules (|A.36p and (IA.13p . respectively for the point particle vertices V,^-^ (j = 
l,---,5) and potential graviton propagator (Hi^^^Hf^^^^), together with relations (|B.1|1 and (1B.3P 
one gets the contributions from each of the diagrams of Fig. [2j 



^Instead, if we prefer to have a Newtonian potential with d- independent pre-factor, Gdm/r'^ ^, we should have 



d-l 

chosen the normalization k^ = 4— i "'"^", o, Gd 

3 d-a r(^^) 
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Starting with the first diagram, one has 



(2) \-'-'K(jaV-^2;-'-'qMi'V'^l// '^(2) 

dt "rf^s^ ^ |„ ,JC"""''/Tm«>-s (^1 • ^2) • (3.6) 



-4 r (^) Gdmim2 

^^ |xi(t)-X2(t)|'^-3 

For the second diagram, there are two point particle vertices, namely Vlt'^ and V,'^ in (|A.36p . scaling 

as L^l'^v^ (see Table [2|). Also, there are two mirror contributions resulting from inter-changing the 
two particles, (1 <-> 2). Taking this into account we get, 

Fig.[2b = V^(g4)(^k„,(:c0)i/q^,(x?))y(^;^ + (l-2) 

'' d^2 ,^ |x,(t)-X2(t)|^-3 (vi + V2) . (3.7) 

To compute the third diagram we introduce a symmetrization factor of ^ to account for the sym- 
metry under the exchange of the two graviton lines, yielding 

Fig.Et = ^l/(if(^k.,(x^)i^q,.(x?))y(^^'''(Fq,,^(x?)/7p,„(x^0))F(^;' + (l-2) 

2 ^^ ; |xi(t)-x2(t)p('^-3)- ^^-^^ 

The fourth diagram requires the Feynman rule for the 3-point correlation function for the interaction 
between 3-potential gravitons, (/fkoo-^qoo-^poo)- This is computed in Appendix IA.5I and given by 
(1A.34|) . The symmetrization factor of ^ is also needed. Using Eqs. (jB.ip and ()B.3p and the integrals 
(1531) we get 

Fig. [Hi = iy(^f(i7k„,(x?)i/q^,(x°)ifp^„(x^20))yJ^y(^;' + (i-2) 




To get the contribution from diagram of Fig. [2^ we need the v^ correction to the potential graviton 
propagator, {H]^^^H^^^)^, given in (|A.15p. We get 

Fig.[2}3 = V^,^{H^^,{xl)H^^AA)Uyi;^ (3.10) 

dtt'^'^'W . .^'"'"l.-. fvi-V2-(d-3)-^"^-"^^^^"^-"^^- 



d-2 ^^ |xi(t)-X2(t)|^-3 V ^ - V ^ |X1-X2|2 

where X12 = xi — X2. To evaluate the result of acting with the Bq operator (appearing in 
(-ffkafl-f^q^^)®) we used the second relation of (jB.2p together with Si^x^ = (^^Vq. 

Summing all contributions (|3.6|) - (|3.1Up and adding also the first relativistic correction to the 
kinetic energy (which is also of order Lv^, as can be seen from the last term in expansion (j2.9p ). 
we finally get 

r _1\^^,,4 o('L:1YSS}\ Glmim2{mi + 1112) 



^"^ ^ '"^-1/ 2 , 2x 3d -5 „Jvi •X12)(V2 •X12) 

vi • V2 — (_a — 6j- 



d-3r(^-) Gctmim2 



d-2 ^^ Ixi-xal'^-s 



vr 2 



a — 3 a — 3 l^i — X2|^ 
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Leih is the Einstein-Infeld-Hoffinann correction to the Newtonian potential in a d-dimensional 
background. For d = 4, (j3.1ip reduces to the result originally derived in [32] and latter reproduced 
with the CIEFT in [10] and [20]. 

A final comment is in order. An improved version of the CIEFT introduced in [TO] was recently 
proposed yL9j. This improved version [1^ is probably the most economic computational way to 
get the desired classical observables when the geometry is stationary, by optimizing the original 
EFT construction. Its key novelty is the observation that for stationary geometries one can do a 
Kaluza-Klein dimensional reduction along the time direction. Generically and in short, this intro- 
duces a scalar and vector fields that replace the original tensorial graviton and whose propagators 
are simpler. This technique also emphasizes the classical nature of the CIEFT and eliminates un- 
necessary quantum features (after all we only want tree-level results). This method has been used 
|20j to recover the EIH correction in 4 dimensions. Most of our computations were done prior to 
publication of [19], which is the reason why we did not follow this approach here. After the release 
of the preprint version of our paper, Kol and Smolkin have derived our results using their improved 
CIEFT Ifj. Their results are due to appear in the latest arXiv version of Ref. |20j . 

4 Quadrupole formula in higher dimensions 

Our next step is to find the Lagrangian L^ad describing the coupling between the NR source and 
the radiation gravitons, up to leading order in v. As we will soon show, this occurs at order Lv 2 
and is done in section 14. 1[ The quadrupole formula is then computed in section 14.21 

4.1 Lagrangian for the interaction between particles and radiation 



To find Lyad we take the small v expansion ()2.9p of the point particle NR ensemble in the long 
wavelength regime, i.e., with 5gfj,i, = h^y. In this expansion we just consider terms that couple 
to the graviton h^j^y. These are the terms that can potentially describe the emission of waves out 
to infinity. Finally, we take the multipole expansion (j2.11|) of the radiation gravitons to organize 
systematically the Feynman diagrams in a ladder of powers of v. 

Consider the first such term in ()2.9p . describing the coupling of /iqo to the mass monopole ^ rua, 



E 



„0'^9i 



n^aj dxl^hoo{x) (4.1) 



T. 



m, 



ij d^aY f ^oo(a;°,X) + -5xj, 9^/100 1 (^o^x) + ^^K^^i^idjhool^^o^:)^) + 0{v^] 



d~3 



Using the power counting rules of Tabled] this expression scales as [dx''][Kgrn,][/ioo] ~ \/Lu~2~(1 + 
V + v'^ -!-•••). We are free to work in the center of mass (CM) frame where one has X = 0. Then, 
the second term in ()4.2p vanishes, and in the last term one has 5x* = x* — X' = x*. The first term 
in (|4.2p gives the leading order contribution to the coupling between particles and radiation, 

4i = -yE"^«^oo(:cO,X) ^ Fig.Hi, (4.2) 

and corresponds to the contribution coming from the Feynman diagram of Fig. [3^, that is Fig. [3^ = 

r (a) , — d+1 

ij dt L^^^. The third term in ()4.2p is stored into the \/Lv 2 contribution. 



^And in the way finding an incorrect factor in the previous version of our final formula (|3.1ip . 
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a) b) 




+ (1«^^2) 



c) d) e) 

Figure 3: Feynman diagrams for the coupling between the NR source and the radiation gravitons, up 
to leading order in v where radiation emission occurs. Diagram a) is of order \jIjV^^ ^ diagram b) is 
O (\/Lv^^ j, and c)-e) are order \/Lv^ . Only external radiation gravitons, represented by wavy lines, 
describe the radiated waves. 



The second coupling term in (j2.9p . 

-"^rua dx^Ugvl^hoiix) (4.3) 

= -^ruaj dx^Ugvl (^/ioi(a;°,X) + (5x^5j/ioi|(^.o_x) + ^i'"'^)) > 

scales as \'Lv^~{l + v + ■■■). By definition, the first contribution is proportional to the CM 
momentum, ^ mavl^ = P^m and vanishes in the CM frame, that we choose to work on. Thus, 

4ad = ^ Fig.Eb, (4.4) 

and the last term in (j4.4p is stored into the \fLv~^~ contribution. The vanishing of the \fLv~^^ 
contribution reflects the absence of radiation emission of dipole nature in gravity. 

The third and fourth terms in (j2.9p both scale to leading order as ^/Lv^2~ . This corresponds 
to taking only the leading order term /i^i,(x'^,X) in the multipole expansion of h^^(x). These two 

terms together with the two terms that have been previously stored yield the following \/Lv^~ 
contribution to the interaction Lagrangian between particles and radiation, 



-^ Fig.Et. (4.5) 

As we shall soon realize, the leading terms in Lrad contributing to radiation emission are of order 

, — d+l , — d+1 

\/Lv 2 . Therefore we disregard the other terms in (12. 9|) . all of order higher than yLv 2 . However, 
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other Feynman diagrams contribute at order yLv^2~ , namely the diagrams drawn in Fig. [3]1 and 
Fig. [3^. They contain both internal potential gravitons and external radiation gravitons. 

The Feynman diagram [3]1 is formally the tensorial product of diagram [1^ with diagram [3^ and 
yields the contribution □ 

Kg 

2 



Fig.EH = \v^,^{H^^^{xl)H^^M))Vl'^\\-'^K,{x\^)\+{l^2) 



i / dtL 



(d) 



L 



id) 



d-3r{^) 



K'imim2 



To find the contribution of the diagram shown in Fig. [3^ we need the 3-point correlation 
function for the interaction vertex between two potential gravitons and a radiation graviton, 
(/i(x'',X)//koo(^i)-^qoo(^2))- This is computed in section lA.41 and given by ()A.30p . Using also 
the vertex rule for the coupling V/?f between a particle and a potential graviton we find|f| 



(1) 



Fig.S. = y"f(^(xO,X)ifk„^(x°)i7,,,.(x?))F(^;' 



(1) 

i I dt 



(4.7) 



Kimim2 d- 3 



d-2 



-ik-(xi-X2)_ 



k4 



k^/iuu_ k^r/.,+k,kjU*^ 



I dt ^2 



) Klmim2 (s^Qo ^ d-3 xi^xj^ ^^ 



L 



(e) 
rad 



d-2 i6^^|xi2(t)|'^-3 \^2'^ ' 2 |xi2(t)|2 ^^^ 
d-3 3r(^) K^r 

a=l,2 



i / dtL 



(e) 



rad 



^-2 32 vr^ |xi(t)-X2(t)| 



h'\x\^)-^ Y, m,xj.i^,/i,,(x°,X) 



where x = 4tt- We used the momentum integrals (|B.3p . (|B.4p and Newton's laws, F^ = rria^a and 



F = —VUn, with C/tv being Newton's potential energy (j3.5p . Indeed, they allow to write 



(d - 3)2 r (^) K^mima xfaX- 



12 



d-2 32^^|xi2(t)|'^-3|xi2(t)| 



12 '"^J 



-^(xl-x^)(V^-^^)/i., 



^ ^ TTlaXlxih 



ij ■ 



(4.8) 



a=l,2 



Summing the contributions from all the Feynman diagrams of Fig. [3l i.e., adding (J4.2p - (j4.7p we 

d+l 

2 . 



get the interaction Lagrangian between particles and radiation up to order \fLv' 

Srad = dt Lrad , Lrad = L^d + ^rad + ^rad + ^rad + ^Tad + " " " ' (4-9) 



-w 



-(e) 



Integrating by parts L^J^ and the last term of -^j.^^, we can write ()4.9p as 



^rad 



V- /. 1 2\ d-3 2r^ 

^mjl + -v2^ ^ 2 



Gdmim2 



,00 



d-2 ^^ |xi(t)-X2(t)|'^-3 



(4.10) 



+ ^ m-a (x>^ - x^v^) dihoj - ^ X] "^ax^x^ {dodihoj + (9o<9j/ioi - ^o/ijj - didjhoo) 



''Using the power counting rules of Tables [T] and [2] diagram [3ji indeed contributes as V(i){-H'kfl'q)V(i)Kg/ioo ~ 

^DiagramU? scales as V(^i){hH\,H^)V(^) ~ ('Li/2i;°/(r''/2u'/2)] ('r'^L-^/^^T^) ~ y/Lv^ , where we used the 
rules of Table [H 
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This Lagrangian can be written in a more suggestive form. Introducing the orbital angular mo- 
mentum, 

Lfc = eijkK^i , (4.11) 

one has J2^a yKi^l - ^l^l) dihoj = e'^'^Lkdjhoi, where we used e'-^^enmk = ^''m^n - 5*„(5{ra. To 
interpret the last term in (j4.1ip . use the definition of the quadrupole moment of a gravitational 
source, 

r^- = J^m,xix^ (4.12) 

a 

and the Riemann tensor Rq^qj of the radiated gravitons (this follows from (IA.2p . taking Minkowski 
as the background metric), 



Kq 



^oiOj = ~^ {dodihoj + dodjhoi - dghij - didjhoo) . (4.13) 

The Lagrangian describing the coupling of radiation to NR sources can then be written as, 



L 



rad 



Y^maU + ^K] -Un 



hoo _ ^e^i^L,a,7io. - 1 r^Rill + • • • ' (4-14) 



where Un is defined in ()3.5p . The first term in Lrad represents the coupling of h^^ to the rest mass 
plus Newtonian energy of the system. It results from the fact that in General Relativity, not only 
the rest mass but also the kinetic and gravitational potential energy contribute to the effective 

/ d — 3 

mass of the system. This mass correction or monopole term is conserved at order \'Lv~2~. Thus, 
it does not contribute to radiation emission, which is a dissipative effect. The dipole contribution 
also vanishes as observed in (14. 4p . 

The second term in Lrad, represents the coupling of the orbital angular momentum of the 
system to the radiation gravitons /iqj, and is of order \/Zf^~. At this order, the orbital angular 
momentum is conserved and thus this term does not contribute also to radiation emission. 

The leading order term describing radiation emission is the third term in Lrad, ~ 2^*''^Oi07' which 

, — d+l 

is of order \/Lv 2 . Note that this term relevant for the emission process arises from diagram l3p, 
which has 3-point graviton self-interaction. This clearly shows that the quadrupole gravitational 
emission is sourced by the non-linearities of General Relativity. Ultimately, this is a consequence of 
the fact that graviton energies also source the gravitational field in General Relativity. In the next 
subsection we use this radiation term to find the emitted radiation power. As in the d = 4 case 
|10] . the Lagrangian (j4.14p is manifestly gauge invariant under infinitesimal gauge transformations, 
showing the power counting scheme in its full glory. 

4.2 Quadrupole formula 

We now want to compute the quadrupole formula that gives the total energy per unit time radiated 
through gravitational waves by a NR source. The total number of gravitons emitted by the slowly 
motion system over the period T is given by 



/ 



■^^dys"!'"'^"-!^-'- c-'^' 
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Here, dT is the differential graviton emission rate, and SqfIxoI is the effective action that follows 
from integrating out the radiation gravitons through the functional integral (j2.13p . The imaginary 
part of this action, IuiSqf, gives the quantity we are interested in and is given by the imaginary 
part of the self-energy Feynman diagram represented in Fig. HI Only the last term in (j4.14p 
contributes to radiation emission since this is the only term with an imaginary part, to leading 
order. The self-energy diagram then gives 

Fig.Sl= -^jdx\dxlr^{x\)I^^{xl) (i?SJ,(^?,X)i?W^(xO,X)). (4.16) 




Figure 4: Self-energy diagram whose imaginary contribution gives the quadrupole formula for the radiated 
power in the form of gravitational waves by a NR system. The double solid lines represent the NR particles. 
In the low energy CIEFT describing the coupling of the particles to the radiated gravitons, lengthscales 
smaller than the orbital distance cannot be resolved. This is pictorially emphasized in this diagram by 
drawing the particle lines close to each other. This diagram is to be understood as a standard self-energy 
diagram: at a certain point in the past the system radiates a graviton (described by (|4.14p ) which is then 
absorbed back in the future. The amplitude for this process is complex. Its imaginary part encodes the 
information on the radiated energy through gravitational waves that propagates to the asymptotic region. 
Again, this is a tree-level diagram since the loop is closed by "heavy" particle worldlines whose momenta 
are not integrated. 



To find the two-point function of the linearized Riemann tensor we take the definition ()4.13p of 

) 

OiOj 



iJgjQ- and the radiation graviton propagator (jA.Sp . yielding 



(41(^?,X)<L(^2,X)) 



ik^ d{d - 3) 



1 .. . ... 1 



2 



y'Jik^jrn ~r Oi-rjiOj]^) (^ij^kva 



(4.17) 



^o^OjV-i,-;-ofcOmV-^.,"// - fc2 + ie4(d + l)(d-2) 

where we used that from the on-shell condition for the radiation gravitons, /cq — k^ = 0, follows 

that kjkj = -^[^5ij and kjkjkfckm = -^p^\5ij^km, + ^ik5jm + ^im^jk)- Note that the two-point 
function in (14.170 is proportional to the projection operator onto symmetric and traceless two- 
index spatial tensors. So when it acts on the quadrupole moment in (j4.16p . only the traceless part 
of the quadrupole moment P^ will effectively contribute. This traceless component is given by 

Introducing its Fourier transform, 

Q*i(A;o) = /"dx°Q^^(x°)e-^'="^° , (4.19) 

and using (j4.17p . (|4.16p becomes 



F.g.il=-a , "^j;',. ik^.i^h:mM\ (4.20) 



t^l d{d - 3) /• d'^k {ky 
32 (d+l)(d-2) J (2^)"! k^ + ie 
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with \Qij\ = Q^-'Qlj- To find the imaginary part of ()4.20p . that gives IthSqf, one uses the 
well-known relation for the principal value of a function, 



^ (^0 -^') = ^\^ (^° -\^\)^ (^° + Ik|) > (4.21) 



where P stands for the the principle value. This relation extracts the imaginary contribution of 
diagram [H 

Here, the integrand gives the differential graviton emission rate along the evolution of the system. 
The differential power radiated is then this emission rate times the energy |k| of the graviton. And 
the total power radiated during the period T ^ cxd is then 

72 



that is 



dEd^E ^^^ = —l,nSQF[xal (4.23) 



"^ ^^^"'^ ^ ^'"^|knQ.,(|k|)|2. (4.24) 



32r {d+l){d-2) J (27r) 
Exploring the spherical symmetry of the problem one has J d k = ^d-2 / (iA;|k| , with il^ 



2 



27r" 



d-1 



^ 



2 



being the area of the unit [d — 2)-sphere. Moreover, the on-shell condition for the graviton 



momentum allows one to make the replacement |k| = —k^. Using k^ = w, and the normalization 
(|3.4p for the gravitational coupling, the energy radiated per unit of frequency finally reads 

^-^' id-2)id + l)T[^] " "^^^^"^'- ^'-''^ 

This is the celebrated quadrupole formula, valid for any spacetime dimension d > 3. For d = 4 
it was originally computed by Einstein and reproduced using the EFT formalism in [10] and |20j . 
For even d, our formula agrees with the expression first computed in [27] using a standard Green's 
function formalism. Using the CIEFT approach we have shown, for the first time, that the formula 
is also valid for odd d-dimensional spacetimes. This was an open problem up to now [27] and this 
is thus one of our main results. 

Relation with previous work 

The reason for the difficulties to derive the quadrupole formula in the traditional approach 
is the following [271 [28]. Wave propagation of massless fields is considerably different depending 
on whether the number of spacetime dimensions d is even or odd. This follows from the fact that 
Green's function has support only on the light-cone in the former case, but also inside of it in the lat- 
ter background. Indeed, in even d the Green's function is roughly [27j G^'^^ ~ dr [r~^(^(i — r) 

" 0{t-r) 



(where r stands for radial distance and t for time), while in odd d one has G^'^^ ~ dr 
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VP 



l^fJLV 



(where 0(x) is the step function). As a consequence, wave tails or wakes develop in odd dimensions 
PH] . Now, traditional methods to compute the quadrupole formula require solving the linearized 
Einstein field equations in the harmonic gauge, D/i^i, = IGttGciS^^, where D is the d-dimensional 
Laplacian, and Sf^u = Tf^u — ^fr2 Vi-iiy ^ encodes the information on the energy-momentum tensor of 
the NR source. The general solution of this inhomogeneous differential equation can be expressed 
with the help of the retarded Green's function G'''^*(t — t',x — x'), as 

(t,x) = IGnGd I dt' f d'^-^x.'Sf,uit',x.')G'^\t -t',x- x') , (4.26) 

plus any solution to the homogeneous equation. For even d dimensions, we have a trivial integral 
over a delta function. However, in odd d, the structure of the Green's function displayed above 
makes the integral in (j4.26p highly non-trivial, and so far this integration has not performed. 

In a forthcoming publication [3D] we will show that the quadrupole formula (|4.25|) can also 
be reproduced using a more traditional Green's function formalism. The key observation will be 
that the Green's function in momentum space has the same functional form for both even and odd 
spacetimes. So, carrying the calculations in momentum space instead of in the position space (as 
is standard) allows to re-derive (j4.25p . We will further discuss in detail the differences in even and 
odd d spacetimes, and explore ()4.25p . 
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Appendices 

A Feynman rules for propagators and vertices 

To compute the Feynman diagram contributions of Figs. [H [2] and [3] we need Feynman rules for the 
graviton propagators, graviton vertices and point particle vertices. In this appendix we give these 
rules and details of their computation. The derivation is self-contained. We adopt the background 
field method introduced by Dewitt [53] in 1967 and fully developed by t'Hooft and Veltman [3^ . 

We start by considering metricperturbations Sg^^, around the unperturbed background gjiJ , up 
to the order we will be interestedjj 

9^.u = gj^J + Kg6g^u , g^'' = g^^^-K,6g^'' + Kl5g'^Jg''' + 0i6g^), 

^We use the (+, —,■■■,—) signature so ^/g — ^/zEdotg^ for odd and even d, respectively. 
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V~9 



m 



^iTriHS^+t.gSg'^,)] 



1 + 



■6g 



^Qap^g' 



afS 



-M 



+ &iSg' 



(A.l) 



where we used g^'^g-yu = S%, Sg = g^Q\5g^u, detM = e'^i''^*^!, and the Taylor expansion for In(H-x) 
and e^. 

The metric perturbation (jA.ip naturally induces perturbations on the affine connections, Rie- 
mann and Ricci tensors. The perturbations introduced in these tensors can be found in [34j and 
a nice reference to find the details of their computation is [39j. Here we present the first order 
perturbation in the Riemann tensor, because we need this quantity to get (|4.13p . 

'-'^Kap = Y {"^cN.Sg^p - Vf^VJg^^^ - VaV^Sg,^ + V^gV^^^.^ +^°^R';^p Sgl +(°)i?X^„ 6g^^) . 

(A.2) 



It also induces perturbations on the Ricci scalar. Expanding the Einstein-Hilbert Lagrangian (j2.7p 
in powers of the gravitational field 6g^,y one finds that L[g] = L^^> + L^^> + L'^^ + 0{5g^) with: 
L^^> = 2w-^'^ / K? , L^^' = (due to the unperturbed Einstein's equations), and 



L(2) 



/5(o) 



2 \2 



Sg' - dg^^Jg''^ R^'^ + {26g'^Jg'''' - 6g Sg^^^ i^W - \ {V^5g) {V^5g) 



1 



+(V^55) {-^u^gn + o (Va%.) (V"^^'"^) - {V^Sgn {"^ M.) 



(A.3) 



after partial integrations (PI) and moduli total divergencies (MTD) in the action integral. 

At this point we have to choose the gauge. The background field method adopts the viewpoint 
in which the background unperturbed metric gf!!^ is left invariant under difeomorphism transforma- 
tions, x^ ^ x^ — (,^ (where S,^ infinitesimal vector field), while the metric perturbations transform 
as dg^i, —>■ Sg^iy + V^^i, + Vj^^^ (where V^ is the covariant derivative w.r.t. gliJ). Hence its name. 



Choosing to work in the harmonic gauge, Va^g°'^ 
Lagrangian and associated ghost Lagrangian [34j . 

1. 



l^fi^g = 0, one then has the gauge fixing 



-^GF = y/g{ff) I VafJs-"^ - -V^6g I Va6g 



^a^l 



1 



V^'55 , 



i^ghct = Vvri*^ (V^V^T?^, - 40J 



r? 



(A.4) 



where 77^ is the fermionic ghost field. The ghost contribution would be important only if we were 
interested in computing quantum loop corrections. Since we only want tree-level results we will 
make no further reference to it. 

Defining L^g2 = L'^^ -|- Lgf one then has 



L, 



6g^ 



\ Q552 - Sg^^g")! i?^°) + (25<<55°^ - 5g 5gn Rf} + \ {V^5g^,) {V^Sg^ 
- {VaSgn (V^5g^) - i (VaSg) {V^6g) + {V^6g»^) {V,5g'^) 1 , (A.5) 
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where the covariant derivatives V^ are taken w.r.t the background unperturbed metric g^^-.. 

At this point, as motivated in discussion before (12. Sp . we take the unperturbed background to 
be Minkowski spacetime, g\tu = Vfj-u- After PI and MTD we then have 



\{do.5g,u)d'^5g^--\ 



L5g2 = ^ (5aH-) d'^Sa^" - - idcSg) d'^Sg . (A.6) 



A.l Radiation graviton propagator 

So far, our discussion is independent on the lengthscale of the gravitational perturbation. We 
now take the graviton perturbation to be a long wavelength radiation graviton, i.e., we make the 
replacement 5gf^u -^ h^y. Then, again after PI and MTD, ()A.6p reads 

L^,2 = -]^Ki,vf^''h^,, with pf^'^ = i(r/-V^ + r?"V^-r?-V'^)5A9\ (A.7) 

and the associated action is 3^2 = J d'^x-L^2 . 

To find the graviton propagator we have to get the symmetric inverse of the bilinear operator 
jyop^iv ^ 'pj^jg amounts to find the inverse tensor P-y^ap that multiplied by the tensorial part of J)'^^^ 
gives the symmetric identity I-ya^'^ defined in (|A.20p . and to find the inverse of the differential 
operator dxd^. This is simply the Feynman propagator for a massless boson, Dp{x — y). The 
Feynman rule for the radiation graviton propagator is then 

aP 'x-x-v/v/x. |av = {ha(s{x)h^„{y)) 

= Df{x - y) Papf,u , (A.8) 



with 



Pap^lu = ^ ( Va,ir]pi, + r]aur]p^^ - j—^ r]apr]^u ) . (A.9) 

2-d,4-2 



This propagator has power counting scale {hafjh^y) ~ r v , as shown in section [2T6l 

A. 2 Potential graviton propagator. Correction to the potential propagator 

To find the potential graviton propagator we go back to ()A.6p . This time we take the graviton 
perturbation to be a potential graviton, i.e., we do the replacement dg^j^u -^ H^y, 

Lh2 = \ {d^H^,) 9"//^- - ^ {d^H) d^H . (A.IO) 

Introducing the Fourier transform (12.10p of -ff^jy, and using the integral representation of the delta 
function (IB.ip one gets 

Lh2 = -\J^ \k'H^^^H^l - Y^k^-kj -\j^ {-d,H^^^d'H^_l + \doH^d'H_^ , (A.ll) 
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where i^k = H^ , and the associated action is 3^2 = J dx^L^2. The terms in the last curved 
brackets are suppressed relative to the square brackets terms by a power of v^. They are treated 
perturbatively, as operator insertions, in correlation functions (see discussion below). So, taking 
by now only the leading contribution of (jA.lip . it can be written as 

Lh2 = -^I^H^^^^V'^/^^H^^^, with P^^^^^ = i (r?- V^ + r,- V^ - r?- V^) k^. (A.12) 

To find the potential graviton propagator we have to get the symmetric inverse of the bilinear 
operator T)'^^^ . The inverse of its tensorial part is PaPfiv, and the inverse of its momentum part 
is simply k^^. The Feynman rule for the potential graviton propagator is then 

aP ^v = {H^^y)H^^M')) 

= -(27r)'^-i<5(k + q)j^<5(xO-x^O)Pa/3M- (A.13) 

with PaPfiu defined in ()A.9P . As shown in section 12.61 the power counting scale for this propagator 
is (i^k ,H„)^r'^v. 

To compute diagram[2^ in the main body of the text we need the next-to-leading order correction 
to the potential graviton propagator. This accounts for the contribution of the curved brackets 
terms in (jA.lip that we neglected to get the leading propagator (jA.lSp . The quickest way to get 
this propagator correction is to note that the graviton propagator is proportional to t^ = tt-tji- 

For a potential graviton, and in the small v limit, one has jt§- = w <C 1 which allows the Taylor 
expansion, 

1 



A;2-k2 



-j^(l + § + 0(.^)). (A.14) 



Since A;^ <-> Bq the potential graviton propagator correction contributing as r'^v'^ corresponds to 
insert the operator 9g (which we denote with the subscript ®) in the correlation function (jA.lSh . 
yielding 



ap - -^ - Hv = (//k.,(x°)i/q,.(x^°))® 

= -{27,f-H{\. + c,f-^5{x^-x'^)P^a.,. (A.15) 



A. 3 3-radiation graviton vertex 



To get the 3-radiation graviton vertex using the background field method we have to go back to 
(jA.SP and expand the background field g^^s^ up to first order, 

afv = Vfiiy + i^g^fiu , (A. 16) 

where the background metric r/^,y is taken to be Minkowski spacetime and hfj^i, represent the new 
small perturbations around it. Of course, ()A.16P also induces perturbations on the metric deter- 
minant y/g((^, affine connections, Ricci tensor RiiJ and Ricci scalar R^^' that appear in (IA.16p . 
The result of this expansion is the Lagrangian for the 3-graviton interaction, 

S^,,2 = I d'xL-,j,2 , Lj,,2 ix) = -^ h^''{x)Tl:l (x) , (A.17) 
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with 



-^d^d^ - - rif,^d^ 






+2da 






h"" d'^hau - h^^d'^haij, + h^^d'^h 



(A.18) 



The Feynman rule for the 3-radiation graviton vertex fohows from its definition in momentum 
space, 



(%2)aLa =i{'^T^f / d'^x / (i''xi(i'^a;2d'^X3e*('=^i+P^2+9X3) 



5h'^P{xi) 6h^-{x2) <5V(x3) 



with 



dhapix) 



5'^{x-Xi)Iai3^iu, k + p + q = 0, 



6ht^^{xi) 

and lajS^u being the symmetric identity tensor 

1 



Ia|3^lv = 2 iVafiTliSu + VauVfitJ.) ■ 



(A.19) 



(A.20) 



We assumed that ah momenta k,p,q is incoming and thus conservation of momentum in this 
convention reads as displayed in the end of (jA.lQp . After a long but straightforward computation 
that also makes use of the integral representation (jB.ip for the delta function we arrive at the 
desired Feynman rule for the 3-radiation graviton vertex, 






+2gAg5 



fc^F + {k + qY{k + qf + qf'q" - ^^"^ q^ 



-^af3-ya „ Vaf^V^c 



r\5 



,I^\rr+I 



XS jnu 



TtJ.S jv\ 



jvX 



J^'K 



tXh 






-2q' 



lafsxs (l'\Jk + qr + ^'%.(fc + qY) - Va<5 (I'^aak'' + / 



jSii 



afS 



k" 



W (l'\pIW + 4/3^^'%. + V'^q'qs ila(sxpl'\^ + I,.Xpl'\ 



+ (A;2 + {k + qf 



I^^aB^W + I^" apl^as'' - oV'"" la^o " TT ^"/^^T'x 



k'l^^pI^\, + {k + qfl^,,I^lp)] 



(A.21) 



This Feynman rule is independent of the dimensionality of the spacetime (contrary to the rule 
for the graviton propagator and potential graviton vertices). In particular, ()A.2ip is the same as 
the Feynman rule first obtained in 4-dimensions by \i7\ 138] (for the details see [39]). The 3-point 
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correlation function for the 3-radiation graviton interaction, {h^i,hoif3h^a), can be obtained trivially 
from the tensorial product of this vertex rule with the radiation graviton propagators (jA.Sp , 

">x.--^' YCT = {ha/3ixi)hf,uix2)h-yaix3)) (A. 22) 



|J,V '^ 



{Kp[x^)h^^{x\)) iy-^u^r^"^" {K,[x2)hi,,[x\)){h,„{x3)h^,{x\)) . 



We do not need (IA.21J) or ()A.22h (we just quote them for completeness), but we will make use of 
(JA.ISP in the next subsections. 

A. 4 3— radiation-potential graviton vertex and correlation function 

The action describing the 3~radiation-potential graviton interaction is 

S-hm = jd'^xLj.H, , L-^H^{x) = -^ h^%x)Tl!^{x) , (A.23) 

where T/J^ (x) is obtained from Tl^ij{x) by doing the replacement h ^ H vn ()A.18|1 . 

We proceed using the Fourier transform ()2.10p and keeping only the leading order terms in the 
velocity. This yields 

<(x)=//e^('^+'^)-T,f(xO), 

+ ^ [\^uK + q«q^ + 2k(„q^) + ^r/^^ (k^ + q2 + 2k • q) j [HfH^^^ - ^ H^H^j 

+ {kiku + q^qn + 2k(iq„)) Hi^H'^^ + (kjk^ + q^q^ + 2k(jq„)) Hi^H'^^ 

-2 (k,k, + q,q, + 2k(,q,)) (^Hl^Hl^ - H^H^^,^ - \ r?^, [H^H^^ - /^k^^)) 

-2 (q,q, + kiq„) H^H^^^ - 2 (q,q, + k^q.) i/^^q.. 

- (k2 + q2 + 2k • q) {H^^.H^^ - H^H^^^) + higher order v terms. (A.24) 

Notice that in this expression we use the Latin letters n, v when only spatial /i, v make a leading 
order contribution, and we use parenthesis () around the indices to denote symmetrization, kucij\ = 
(kjqj +qikj)/2. 

Proceeding, we do the multipole expansion ()2.11|) of h^u[x) around the system's center of mass 
X, and keep only the leading order term in the velocity v, h^y{x^ ,'K). So, to leading order /i^j/(x) 
only depends on x^ but not on x. This allows to write ()A.23p in lowest order as 

S-nm = j dx° V(xO, y^)j^j j d'-'^ e^^^^^y- ("f T^f (x°)) . (A.25) 
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The spatial integral can be done using the integral representation of the delta function J e*^^'^-*'' = 
5(k + q). This operation results in q — > — k in ()A.24p and only the first line survives to the 
procedure. We can finally write the leading order action for the 3-graviton interaction between 
potential and radiation gravitons as 



'/l_H"2 



/"(ix°L;i^2(x°,X), with 



L^^.(X^X) = Kgh^' / k^ 



OJ 



+ Kgh 



+ Kgh'^ 



^ko-^-k^o - ^H^^oqH^]^ - -..^ 



2i^i H^\i — H\i,. H—\c 



Hf'H^i^^^ + --ffk^-k 



1 k^ / 1 

-k^i^k ^-k Vii i^r-H^-k - -Hi^H^y, 



(A.26) 



where h^^,h°\h^^ are functions of (x°,X) and Hl^" = i/r(^°)- 

The Feynnian rule for the 3~radiation-potential graviton vertex follows from the definition, 



in 



and use of the functional derivative, 



i dx^ dx\dxl 



5 



5L-,H^ix^,X) 



<5<(x;) <5<(x^) ' 



(A.27) 



J(k + q),5(x°-xO)/, 



i ) ^af3fi.i/ ) 



(A.28) 



with Ia/3^u defined in (|A.20p . Finally, the 3-point correlation function for the interaction between 
two potential gravitons and radiation graviton is (dropping the bars over the graviton) 




|iV 



(/i(xO,X)/7k„,(x?)Fq^,(xO)) 






mg(27r)'^-^5(x^ - x^)5(k + q) 



jj2q2 






1 2 

2 PafS/iiy + PaPt^O VOu + -fa/30i/ VOfi — j _ n (-fa/300 ViJ-u + Val3 PoOfxi/) 



2 {PalBfii riOv + PaPiu VOfi) + , _ „ (??a/3 -PoJAii/ — -fa/30J ??Mi/) 



(A.29) 



'« j "•" 9^*i / ^af3^u ~r k I iafi^j 'Hiy \ ^af3ji/ fJitJ. t q 



^ctj3ij VfjLV 
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In the main text we need only the a, /3, /i, z^ = component which reads 

(/l(xO,X)/7koo(^?)^qoo(^2)) 



(A.30) 



-iKg{2Trf''5{x1 - x^2)S(k + (I 



1 d-3 



k2q2 d-2 



2,00 



■k'h' 



1-2 



k^rjij + kikj W^ 



d+3 



Using the power counting rules of Table [T] we find that {hH]i^„Hf^^J) ~ L ^''^v 2 

A. 5 3-potential graviton vertex and correlation function 

To obtain the Lagrangian for the 3-potential graviton interaction we begin by making the replace- 
ment h^'^ -^ H^'^ in ()A.23p . We then take the Fourier transform (j2.10p of the potential graviton 
H'^'^ and use the integral representation of the delta function (jB.ip to get 



''HH'2 



dx''L 



HH'2 



LHm{A = -^ 



0\n-H'^ lJ^\ 



k Jq Jp 



(2^r^5(k + q + p)i7^^(x")T^^, {x 



(A.31) 



where T^J {x^) is defined in (jX24]) . 

The Feynman rule for the 3-potential graviton vertex follows from the definition, 

5 5 SLhh^x^) 



(yHH^)af^,u,a = i I dx' I dx^dx^^dx^ 



6Hf^{x\) 6H^M) ^K:i4) 



(A.32) 



and use of ()A.28p . The 3-point correlation function for the interaction between three potential 
gravitons is finally 



a|3 



|iV 



--ya = {H^^Ax^i)H^,A4)Hp,A4)) 



(A.33) 



{H^^,{xi)H^^.(x\^)){VsH, 
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±{2^f-'6{xl - xl)5{x\ - xl)6 (k + q + p) j-^-^ 



1 



2P^\^ ( PaPvi/^ Vu)X + -TZr^ P\mP Vf^u ] + J3^ Pap,,u Vya (k^ + q^ + 2k • q) 



d-2 
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d-2 
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where we used (again) the rule that a sum over Latin indices means that only spatial indices 
contribute (to leading order), and the parenthesis () around the indices denote symmetrization. In 
the main text we need only the a, /?, /x, zv, 7, a = component which reads 

{ifkoo(^?)^qoo(^2)^Poo(^3)) (A.34) 

= .,(2vr)^-M(x? - 4)6(4 - 4)S (k + q + p) 1^ ^i$^- 

Using the power counting rules of Table[T]we find that {H]^^ Hq ^Hp ^) ~ L^^'^w^'^r^'^'^. 

A. 6 Point-particle vertices 

The Feynman rules for the non-linear vertex interactions of particles with potential gravitons are 
obtained by taking the small velocity expansion (j2.9p of the point-particle action 5*^^, with 5g^u — > 
H^y. In the expansion (|2.9|) we only display the lowest order terms that contribute to the leading 
order observables of interest. We also take the Fourier transform (j2.10p of the potential graviton 
to keep a well-defined track of the velocity power counting. 

The Feynman rules for the interaction vertex between the particle and a potential graviton 
vertex then follows from the definition of vertex interaction factor, 

and use of ()A.28p . The point-particle vertex Feynman rules are then (following the order of ap- 
pearance in (j2.9p ) 



F« = -^|<|^e"^- 



Tio^imu ■ 



yfj = -^ / < I e^'^- v: ivo.m. + VouV^,) , 



k 



Vi'J = --^ / dxl / e^*^- v[^i {^.,r,,y + rj^yrj,,) 






^1/ 



k 
k 



^2, 



vHtip = ^ J dxl j^ j e^^^+^y^ m^VouVoaVop • (A.36) 

The power counting rules for the vertices are shown in Table [2] (the kinetic contributions in (j2.9p . 
i.e., the two last terms, have power counting rules of L^''^v^ and L^'^u^, respectively). 

B Useful relations for Feynman diagram computations 

In this appendix we present some standard formulas that are useful to compute Feynman diagram 
contributions. 

We start with the integral representation of the delta function 

5(k) = 6''-\k) = J -^^e'"^-^ EE ^ e^k.x ^ ^(^) ^ yT g,,k.x_ (B 1^ 
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and the well known integrals, 

J f{z)5{z -a) = f{a) , J f'{z)6{z - a) = -f\a) , (B.2) 

where f'{z) = -J^. 

The integrals over the momentum are computed using the following relations, 

d"k 1 k.x 1 r(n/2-a) /x^y-'-/^ 



(27r)" (k2)" (47r)"/2 r(Q) 



d'^-ik k,k, .J,., _ rj^ (,_3) 

■ ^ ~ d-1 w 



(27r)^-i k4 8^^ 



5.,-(d-3)^ 



(B.4) 



and 



d'^k 1 2,-2 _r(2-r/)[r(r?-l)]^ 

~ Jo (P ) , -'0 = 



(27r)2^ k2 (k - p)2 "^"^ ' ' ' " (47r)^r(277 - 2) ' 

rf'^k k, _ 1 2N'7-2 

(27r)2^k2(k-pf"2'°P^^P^ ' 

d k kjkj Oij , 2\»?-i , ^ T I 2\»?-2 /o r-\ 

(2vr)2.k2(k-pf ="4(2;FT) °^P^ +P^P^-2(2;FT) °^P^ " ^ '^ 

The integral (jB.SP is needed to evaluate all Feynman diagrams. Relation ()B.4p is needed to compute 
(|4.7p . and relations (|B.5P are needed to compute (j3.9p . 
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